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ABSTRACT 

We perform a cosmological parameter analysis of the 75 square degree CTIO lensing survey in conjunction 
with CMB and Type la supernovae data. For ACDM cosmologies, we find that the amplitude of the power 
spectrum at low redshift is given by erg = 0.811q' X q (95% c.l.), where the error bar includes both statistical 
\^q". and systematic errors. The total of all systematic errors is smaller than the statistical errors, but they do make 

■ up a significant fraction of the error budget. We find that weak lensing improves the constraints on dark energy 

as well. The (constant) dark energy equation of state parameter, w, is measured to be -0.89±£^ (95% c.l.). 
Marginalizing over a constant w slightly changes the estimate of erg to 0.79^0 14 (95% c.l.). We also investigate 
r*i ! variable w cosmologies, but find that the constraints weaken considerably; the next generation surveys are 

D ■ needed to obtain meaningful constraints on the possible time evolution of dark energy. 

Subject headings: gravitational lensing; cosmology; large-scale structure 

1. INTRODUCTION 

Observations of weak gravitational lensing, the coherent distortion in the images of distant galaxies, hav e advanced rapidly in 
• the past four years . The first detections of weak lensing in blank fields were reported only a few yea rs ago JWittman et q/.l2 000: 
^-j- | iKaiser et q/.ll2000t IVan Waerbeke et a/.ll2000t IBacon et a/.ll20061 iRhodes. Refregier. & G roth 200 0]). More recent lensing mea- 
{SJ , surements dHoeks tra et a/.l[2002t [Refregier et a/.l[2002 t IBacon et a/.ll2003t iBrown et a/.ll2003l lHamana et a/.ll2003t Darvis et al.\ 
£NJ ■ 120031 IVan W aerbeke et aZ.ll2005t iHevmans et aZ.ll2005l) have used larger and/or deeper surveys to reduce the statistical errors, 
CD which scale as A -1 / 2 where N is the number of galaxies measured. Better techniques for reducing systematic errors have also 
been developed, resulting in interesting cosmological constraints from lensing surveys. 

Shear correlations measured by lensing surveys determine the projected power spectrum of matter fluctuations in the Universe. 
These fluctuations are believed to have grown due to gravitational instability from the early universe to the present. The growth 
O of fluctations from last scattering, z = 1100, to the present is sensitive to the densities of dark energy and matter via the Hubble 
expansion rate. Further, the measured lensing signal depends on angular diameter distances to the source galaxies. Thus weak 
lensing observables probe the dark energy density through both the growth function and the geometric distance factors. Present 
lensing measurements are sensitive to redshifts < z < 1. 

As the size of weak lensing surveys increases and the statistical errors keep going down, it becomes more important to similarly 
reduce the systematic errors. There are a few different systematic errors which can contaminate a weak lensing signal at the level 
of typical cosmic shear measurements. The largest of these has typically been the corrections of the anisotropic point-spread- 
function (PSF). We have recently developed a principal component analysis approach to interpolating the PSF between the stars 
in the image using information from multiple exposures (Jarvis & Jain 2005). The PSF pattern is found to be a function of only 
a few underlying variables. Therefore, we are able to improve the fits of this pattern by using stars from many exposures. We 
have applied it to the CTIO survey data presented in lJarvis et al.\ d2003l) and have shown that the level of systematic error is well 
below the statistical errors. Indeed the measured B-mode in shear correlations is consistent with zero on all scales measured. In 
this paper we present a new parameter analysis of the CTIO lensing survey. With the reduced systematic error, we are able to 
extract information from the measured signal over nearly two decades in length scale. 

Type la supernovae observations led to the discovery of the accelerated expansion of the universe (Riess et al. 1998; Perlmutter 
et al. 1999). By combining information from the CMB at z = 1100, large-scale structure and SN Ia, interesting constra ints on 
dark energy have been obtained ( Sperge l et al l2003tlBridle et a/.l2003HWeller & Lewisl2003tlAlam et q/.l2004|ISaini et a/.l2004t 
Tegmark et al. 2004; Wang & Tegmark 2004; Seliak et al. 2005; Simon. Verde. & Jimenez 2005; Rapetti et al. 2005; Jassal et al. 
2005). Whether the dark energy density is constant or evolving with cosmic time is one of the most interesting observational 
questions. It is often expressed using the parameterization 

w(a) = w + w a (l - a), (1) 

where w is the dark energy equation of state parameter w = p/p and a = 1/(1 + z) is the expansion scale factor 
JChevallier & Po larski 2001; Linder 2003). The cosmological constant corresponds to wq = — 1, w a = 0. A non-zero value 
of w a corresponds to a time-dependent equation of state. We investigate three priors for w: ACDM (wq = — 1, w a = 0), constant 
w (w a — 0) with —3 < w < 0, and variable w with —8 < w < 8 and —8 < w a < 8. 

In <2] we briefly review our CTIO s urvey data and shape measurement technique, largely deferring to our previous papers 
dJarvis et a/.ll2003l iJarvis & Jainl l2005) for more details. In <J3] we present the results from the new reduction. We use these 
results to constrain cosmological parameters in ^1 an d conclude in |5] 
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2. DATA 

Our CTIO survey data was originally described in detail in lJarvis et al\ ( 120031) . and we refer the reader to that paper for most 
of the details about the data and the analysis. Here, we present a brief summary of the data set, and point out two significant 
changes in the analysis: our PSF interpolation and the dilution correction. 

The data were taken at Cerro Tololo Interamerican Observatory (CTIO) in Chile from December, 1996 to July, 2000. We 
observed 12 fields, well separated on the sky, in low extinction but otherwise random locations. Each field is approximately 2.5 
degrees on a side, giving us a total of 75 square degrees. The 50% completeness level occurs near R = 23.5 for each field, 
although it varies somewhat between the fields. We impose limits of 19 < R < 23, which gives about 2 million galaxies to use 
for our lensing statistics. 

The shape measurements of the galaxies follow the techniques of Bernstein & Jarvis (2002). The galaxy shapes are measured 
using an elliptical Gaussian weight function which is matched to have the same ellipticity as the galaxy. (That is, we use a circular 
Gaussian in the sheared coordinate system in which the galaxy is round.) The observed ellipticity is then calculated as: 

e = e x + ie 2 = - ^ V ^ W ^ X ' y ^ x + (2 ) 
/ I(x, y)W(x, y)(x 2 +y 2 )dxdy 

where / is the intensity, W is our Gaussian weight, x and y are measured from the (weighted) centroid of the image, and the 
boldface indicates that e is a complex quantity. 

We correct for the effects of the point spread function (PSF) in two steps. First, we correct for the effect of the shape of the 
PSF by reconvolving the observed images with a spatially varying convolution kernel which is designed to make the stars round. 
The galaxies in the convolved images are thus no longer affected by the shape of the point spread function (PSF), but are still 
affected by the size of the PSF. 

Since the convolution kernel is only measured where we have an observed star, and the PSF is far from uniform across each 
image, we need to interpolate between the stars. In our previous analysis, we used a separate interpolation for every image. For 
this analysis, however, we use our new principal component algorithm which u ses the information from all of the images at once. 
This new method, which gives a better fit, is described in Jarvis & Jain (2005). 

The effect of the size of the PSF is called dilution. A perfectly round PSF blurs the images of galaxies, which reduces the 
observed ellipticity. For a purely Gaussian PSF and Gaussian galaxies, the measured ellipticities are reduced by a factor R: 

R = 1 - ^ (3) 

However, galaxies are certainly not Gaussian, and stars are only approximately Gaussian. In our previous analysis we used 
a formula intended to account fot the first order corrections due to the kurtosis of the galaxies an d the PSF. Unfortunately, our 
formula was incorrect with respect to the PSF kurtosis, as pointed out by Hirata & Seliak (2003). They give a more accurate 
correction scheme which accounts for both kurtoses correctly to first order - we implement this scheme in the analysis presented 
below. Their formula is quite complicated to write, but is easy to implement, so we defer to their Appendix B for the relevant 
equations. 

Finally, we estimate the shear, 7, from an ensemble of ellipticities using the formula: 

7=^^ (4) 

where the responsivity, 1Z, describes how our mean ellipticity changes in the presence of an applied distortion. It generalizes the 
formula given in equation|3]to non-Gaussian shapes. The factor of 2 in the denom inator above con verts the ellipticity to shear, 
and w is our weight function. We use the "easy" weight function given in Bernstein & Jarvis (2002) (Equation 5.36): 

w = [e 2 + (l.5a v ) 2 }-V 2 (5) 

where a v is the shape uncertainty in th e sheared coordinates where the galaxy is circular. The corresponding responsivity, 1Z, is 
also given in lBernstein & JarvTsI (12002) (Equation 5.33). 

3. WEAK LENSING STATISTICS 

To describe the two-point statistics of our s hear field, we use the aperture mass statistic dSchneider et a/119 98: Critt enden et al.\ 
l200llSchiieider et a/.l200llPen et alJ200l . 

The predictions from theory come in the form of the convergence power spectrum: 

where x is the radial comoving distance, xh is the horizon distance, r(x) is the comoving angular distance, a is the scale size of 
the universe, P^d is the three-dimensional power spectrum of the matter fluctuations, and 

g(x)= r d X ' P (x') r(X '~ X) (7) 

J x r (x') 

where p is the normalized (to give unit integral over x) redshift distribution of source galaxies. These predictions are not reliable 
for high values of k (k > 10/iMpc" 1 ) ( Sm ith et at . 2003) due to difficulties in predicting the non-linear growth. 
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The observed second moments are completely described using the two-point correlation functions: 

£+(0) = <7(r)7*(r + 0)) 

e- (0) + i^x (0) = ( 7 (r)7(r + 9)e- 4lQ ) (8) 
where 8 — 6e la is the separation between pairs of galaxies, and treating the positions on the sky as complex values. The 
correlation functions are not measured at all on scales larger than the size of the survey fields (for this survey, at 9 > 200'), which 
correspond to low k values. 

Using the information from the correlation function to obtain the power spectrum, or vice versa, requires an extrapolation away 
from the k values which are well measured or well predicted. The aperture mass statistic is in some sense a compromise between 
these two statistics, since it can be calculated from either the power spectrum or the correlation function using only the range of 
k values which are either predicted or measured, respectively: 



(mL)(r) = ±- I £dep K (£)w(m) 2 



u.«.r_(!j + ^)T 



(9) 



ap/V , 

_i rede 

~2j & 

where we use the form suggested bv lCrittenden et a/1 fc002l) for which we have 

W( V ) = ^e-^ 

+ W " 128 

T_(z) = ^e-* 2 /4 (10 ) 

The aperture mass therefore has both good predictions from theory and accurate measurements from the data. 

While the integral in Equation|9]is technically from to infinity, the T functions drop off very quickly, so that the effective 
upper limit is really around 5R. Our fields are 200' along the long diagonal, so we can measure the aperture mass up to R = 40'. 
The lower limit, due to difficulties of measuring the correlation function on very small scales, is about R = V. 

The function W{£R) is a narrow function around £R = y/2, so this corresponds to a range in £ of approximately 120 < £ < 
5000. 

The other big advantage to the aperture mass statistic is that it provides a natural check for systematics. Weak lensing should 
only produce a curl-free E-mode component, so any _B-mode observed in the shear field represents a systematic error of some 
sort. The aperture mass statistic has a £?-mode counterpart which can be likewise calculated from the correlation functions as: 



{Ml){R) = \J 6 -§ 



(11) 



Figure ^ shows the results for our reanalysis. The blue points are the E-mode signal, and the red points are the S-mode 
contamination. Points separated by at least one other point are very nearly uncorrelated. The black curve is the best fit flat 
ACDM model found below. 

The B-mode is seen to be consistent with zero at all scales, which was not the case for our previous analysis. Further, in Jarvis 
& Jain (2004) we show that the measured ellipticity correlation function of stars, which is another measure of systematic errors, 
is one to two orders of magnitude smaller than the expected lensing signal at all scales. Therefore, we can now confidently use 
all of the aperture mass values from 1' to 40' for our constraints on cosmology. 

In addition to the aperture mass, we also measure the variance of the mean shear in circular apertures: 

AJx(£R) 2 



(\l\) 2 (R) = ^ J £d£P K (£)- 



(£Rf 



2R ede 



ri^ s +{r) (12) 

where 

' ' ' (13) 



S+(x) = — (4arccos(a;/2) - x\f A - x 2 ^j 



Figure [2 shows the results for the shear variance statistic. There is no E/B decomposition for this statistic 1 . We assume 
implicitly that there is no _B-mode contamination in the shear variance measurements, which seems reasonable given the low 
£>-mode seen for the aperture mass. 

The usefulness of this statistic is that it is able to probe the power spectrum at somewhat smaller £ values than the aperture 
mass statistic. The upper limit in the integral in Equation^]is only 2R, so we can calculate the shear variance up to R = 100' 
with our data. This probes the power spectrum down to £ of about 70. This leads to a total dynamic range for both statistics of 
almost 2 orders of magnitude. The shear variance below R = 50' is degenerate with the aperture mass, so for our constraints, we 
only use the shear variance at large values of R where it is providing extra information. 

In Figure ^ we a ls° show the overall best fit ACDM model (see §3). This fit has a x 2 of 7.7, for effectively 6 degrees of 
freedom, yielding a reduced x 2 of 1.28. 



1 Technically, there is i Schneider et al. 2002), but it requires extrapolation of the correlation functions past where they are measured. 
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Fig. 1 . — The aperture mass (left) and shear variance (right) measurements for our CTIO survey data. For the left plot, the blue (upper) points 
are the iJ-mode measurements, (M% P )(R) and the red (lower) points are the B-mode contamination, (M^){R). For the right plot, the blue 
points are the total shear variance, including any potential B-mode contamination. In both cases, the black curve is the best fit flat ACDM 
model found in i|4] 



4. ANALYSIS 
4.1. Dark Energy Constraints from Weak Lensing 

Our lensing measurement constrains the shear power spectrum, which is a weighted projection of the mass power spectrum. 
The constraints on dark energy arise from two sources. The first is the angular diameter distances to the lens, to the source, and 
between the lens and the source that enter into the projection. The second is the amplitude of the power spectrum. The dark 
energy component alters the expansion rate of the universe at redshifts below about 2 (at least if its evolution is not too different 
from a cosmological constant). This in turns affects the growth of structure. Since the CMB fixes the amplitude of the power 
spectru m at z = 1100, the lensing measurement of the amplitude at low redshift measures the growth function. See lHu & Jainl 
(2004) for a detailed discussion. 

We assume that massive neutrinos make a negligible contribution to the matter density, that the primordial power spectrum 
index has no running and that the universe is spatially flat. The shape of the mass power spectrum is then specified by the baryon 
density tlhh 2 , the matter density il m h 2 , and the primordial spectrum. Following the WMAP convention, we use the scalar 
amplitude A s and spectral index n such that the shape of the primordial power spectrum is A s (k/kop n ~ 1 ' , where fco = 0.05 
Mpc -1 is the normalization scale. The current uncertainties in these parameters are at the 10% level or better. Thus the power 
spectrum as a function of k in Mpc^ 1 (not h Mpc -1 ) in the matter dominated regime can be considered as largely known. 

The amplitude of the power spectrum at a given redshift depends on the initial normalization A s and the "growth function" G 
defined by 



P(k,z) 



1 G(z) n2 



l + z G 



P(k,0), (14) 



where Go = G(z = 0) and we assume that all relevant scales are sufficiently below the maximal sound horizon of the dark 
energy. 

The normaliza tion of the linear power spectrum today is conventionally given at a scale of r — 8/i _1 Mpc and can be approxi- 
mated as llHu & .Tainl2004l) 



^8 



fn h h 2 \- Q - 33 fn m h 2 \ 0563 , n 1( _ Nf „_ n/2 / h \ 0693 Go 



0.97 V 0.024/ V 0.14 7 v ' ' V0.72/ 0.76' 



x (3.123/i) (n - 1)/2 hr^ w^, (15) 



Thus a measurement of as, in conjunction with constraints on the other parameters in the above equation from the CMB, con- 
strains a combination of dark energy parameters that affect Go- Note that while the equation above illustrates how the dark 
energy parameters are linked to others, we do not actually use it in our parameter analysis. Instead we use the projection integral 
of Equation|6]which includes the full range of redshifts probed by our survey. The peak contribution to the lensing correlations is 
from z ~ 0.3, though the maximum sensitivity to w is at z ~ 0.4 as discussed below. Combining our measurement with others 
such as Type la supernovae, which are sensitive to different redshifts, enables a probe of the time dependence of dark energy 
parameters. With photometric redshifts, this could be done with the lensing data alone using the auto and cross-spectra in redshift 
bins. 
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The dark energy modifies the expansion of the universe according to the equation (for a flat universe)(e.g. Lind er & Jenkins! 
12001 : 



H\a)=H* 



(16) 



where the dark energy density is S7dc( a ) = 8TrGpd c /3H(a) 2 , its equation of state is w(a) = Pdc/ Pdc, an d we indicate the present 
time values, fl m (a = 1) and ftd c { a = 1) as simply il m and fl de . 
Distances are then given by 

, , f 1 dale 

X{a) = / TT( 7 72 < 17 > 

J a H(a')a' 

and the growth function G depends only on the dark energy throught the equation: 



d 2 G 



- - |w(a)fi de (a) 



dC 3 

[1 - w(a)]n dc (a)G = . (18) 



din a 2 

When we use a constant w parameterization, it is equivalent to a measuremen t of w at the redshift for which the errors in the 
constant and time-depen dent p iece (in a Taylor expansion) are uncorrelated. See lHu & .Tainl ll2004t) for a discussion of this pivot 
redshift, which we find ( fl4.2.3t to be about 0.4 for our survey (when combined with the CMB and SN priors as described below). 

4.2. Joint Constraints 

We use the results from WMAP ( Soergel et al. 2003) as priors for our analysis. Specifically, we use the Monte Carlo Markov 
Chain (MCMC) calculated by Ver de et aZ.l li2003) 2 . We choose to use the less restrictive prior of it; > —3 rather than w > — 1 in 
order to be as conservative as possible. However, we do have a hard prior that ilk = 0. While WMAP has constrained this to be 
0.02 ± 0.02 for pure ACDM, and there is theoretical bias in believing it is exactly 0, it is worth emphasising that our dark energy 
constraints would be weakened if this prior were relaxed. 

Each step in the Markov chain contains a value for each of the following parameters: cjb = Qbh 2 is the density of baryons; 
uj c = Q, c h 2 is the density of cold dark matter; 6>a is the angular scale of the acoustic peaks; n is the spectral slope of the scalar 
primordial density power spectrum; Z — exp(— 2r) is related to the optical depth (r) of the last scattering surface; A s is the overal 
amplitude of the scalar primordial power spectrum; h = Ho / '100 is the Hubble constant; w is the equation of state parameter for 
the dark energy. 

The parameters 6a and Z are not directly relevant to the aperture mass statistic which we have measured, so we marginalize 
over these two. The others define a cosmology from which we can predict the aperture mass statistic on the scales where we 
have measured it. We present results for selected parameters below, these are obtained by a full marginalization over all the other 
parameters listed above. 

For the linear power spectrum, we use the tran sfer function of | Bardeen et al\ lll986l) . We then estimate the non-linear power 
spectrum using the halo-based model of model of Smith et al. (2003). Their fitting formulae provide a means of estimating the 
quasi-linear and non-linear halo contributions to the power spectrum based on the linear value and the effective spectral index. 
This model agrees with the results of A-body simulations better than the simpler formula of Peacock & Dodds ( 1996). The 
nonlinear correction affects the predict ed variance i n the mass aperture statistic on scales below about 4 arcminutes. 

We also note that the model of S mith et p f] ( 2003 ) th at we u se d oes not include w dir ectly. We correctly take it into account for 
the growth factor and the distances, but Mainini et al. (2003) and Klvpin et al. (2003) show that dark energy changes the virial 
density contrast, A c , which results in changes in the power spectrum at high k values (k > 1 /iMpc -1 ). However, the effect is 
smaller than the expected error in the non-linear model even for a w ~ —0.5 cosmology. 

From the predicted power spectrum, we calculate the aperture mass and shear variance statistics using Equations [6][9] and 1 121 
Our data then give a likelihood value for each cosmology, which we com bine with the CMB likelihood from the MCMC. 

We also use the recent supernova measurements of lRiess e t al. (2004) to further constrain the results. For this, we use the % 2 
calculation program of lTonrv et al.\ d2003l) 3 . 

4.2.1. ACDM Models 

With dark energy priors of w = — 1 and w a = 0, the likelihood constraints are non-trivial contours through a five-dimensional 
parameter space. However, most of the gain in constraining power from the addition of the lensing data comes in the quantities 
fi m and as- We show the error contours projected onto the O m — erg plane in Figure|2] The left plot shows the contours starting 
with the CMB data set, and sequentially adding the supernova and lensing data. The right plot shows the contours in the f2 m — <r$ 
plane separately for each of the three data sets to indicate the degree of their complementarity, which is why their combination 
leads to the tight overall constraints. In both plots, the contours correspond to 68% and 95% confidence regions (A^ 2 = 2.30 and 
6.17). The crosses are at the peak likelihood in the projected plane, which is at: (f2 m = 0.26, as — 0.82). 

4.2.2. Constant w Models 

We show the uncertainty contours for the dark energy priors of — 3 < u>o < and w a = in Figure [3] The left plot give the 
projection in the O m — ag plane, which indicates that allowing w to be free does not significantly worsen the constraints on Q m 
and <jg compared to the pure ACDM model. The fide — w plot (right) shows why. While none of the three data sets individually 
have tight constraints in this plane, the combination of all three leads to a fairly tight contour near (and consistent with) w = — 1. 
The peak likelihood models in the two projections are: (f2 m = 0.25, as = 0.79) and (fide = 0.75, w = —0.90). 

2 Available at http://www.physics.upenn.edu~lverde/MAPCHAINS/mcmc.html 

3 Available at http://www.ifa.hawaii.edu~jt/SOFT/snchi.html The code there was modified slightly to read the data of Ri ess et al\ 120041) . 
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Fig. 2. — Contour plots of \ 2 f° r the (fim, fs) plane. The left plot shows the effect of adding the data sets sequentially, starting with the CMB 
constraints, then adding the supernova and lensing data. The right plot shows contours for each of the three data sets separately. In each case 
the contours enclose the 68% and 95% confidence regions. The x is the best fit model. All other parameters are marginalized over as discussed 
in the text. 




Fig. 3. — Contour plots of \ 2 f° r the (fi m , as) plane (left) and the (fide, to) plane (right) for the constant w dark energy models. Both plots 
show the effect of adding the data sets sequentially. In each case the contours enclose the 68% and 95% confidence regions. The black x's are 
the best fit models in each plane. The cyan x in the left plot is the best fit from the ACDM prior (Figure^. 



4.2.3. Variable w Models 

Finally, we consider dark energy priors of —8 < wo < 8 and —8 < w a < 8. It turns out that some of the dark energy models 
in this range have Qde(2 = 1100) w 1, That is, the mass-e nergy of the universe was essentially all dark energy at the epoch 
of rec ombination. This seems to be ruled out by WMAP data (Cald well et a/.ll2003tlCaldwell & Dorarl2004llWang & Te emark 
2004). Therefore, we make the additional prior that £1de(z = 1100) < 0.5. In practice, all the models have Qde(2 = 1100) ~ 
or 1, so this contraint is effectively 17de(^ = 1100) « 0. 

Given this constraint, the primary effect of dark energy on the CMB is through the distance to the last-scattering surface, d^ss- 
Therefore, we approximate the CMB likelihoods by using the WMap constant-w Markov chain mentioned above, modifying the 
dark energy parameters to maintain a constant rfLSS- Specifically, for each line in the Markov chain, we determine c^lss from the 
values of (l m and w; we select w a from —8 < w a < 8; then we determine what wq with this w a and the same fi m maintain the 
given value of g?lss> and we write these values out as a line in a new pseudo-chain. 

The main approximation in this process is that we neglect the difference of the integrated Sachs-Wolfe (ISW) effect between 
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Fig. 4. — Contour plots of \ 2 for the (fi m , as) plane (left) and the (wo,w a ) plane (right) for the variable w dark energy models. Both plots 
show the effect of adding the data sets sequentially. In each case the contours enclose the 68% and 95% confidence regions. The black x's are 
the best fit models in each plane. The cyan x 's in the left plot are the best fits from the ACDM and constant w priors (Figures[2|and[3j. 



the two models. There is some indication 4 that for some of the models that fall within our contours (Figure^ right), the ISW 
may spread the first peak in the power spectrum of the CMB enough to disfavor these models by 1 or 2 sigma. However, we defer 
a more detailed study of this effect to future work. 

The error contours projected onto the fi m — erg plane and the wo — w a plane are shown in Figure^ With this much freedom 
in the models, we do finally lose much of our constraining power. Even with all three data sets, the error contours are still quite 
large. The fl m — as plot (left) shows that allowing w to vary with cosmological time does significantly worsen the constraints, 
allowing much lower values of erg. 

Also, while lensing does significantly reduce the allowed parameter space in the f2 m — as plane, this reduction has only a small 
effect in the wq — w a plane, so that our constraints there are not much better than those for the CMB + SN data alone. The data 
are consistent with a cosmological constant (wq = — 1, w a = 0) at the 95% confidence level, but w a may range as high as 2 5 . 

We can use the direction of the contours in the wo—w a plane to determine the redshift at which we have the strongest constraints 
on the dark energy. If we change variables from {w(z = 0), w a } t o \w{z = 0.4), w„\, the likelihood contour becomes ro ughly 
vertical. This in dicates that our pivot redshift, or "sweet spot" (Huterer & Turnej l2001l IWeller & Albrechnl2Q02t lHvJl2002t 
iHu & .Tainl2004l) . where the constraints of the dark energy are strongest, is at a redshift of about 0.4. (Of course, the banana shape 
makes this impossible to do precisely, so z p j V = 0.4 is just an approximate value.) We constrain w at this redshift, marginalizing 
over w a (and everything else) below. 

4.3. Systematic Errors 

Systematic errors are harder to estimate than statistical errors, since by their nature they contaminate the data in unknown 
ways. There are four systematic errors which we investigate and attempt to estimate: residual anisotropic PSF as estimated by 
the residual S-mode, calibration error, redshift distribution error, and errors in the non-linear prediction. For shorthand, we refer 
to these as B, CAL, Z, and NL respectively. 

When estimating the contributions of these systematic uncertainties to our error budget, we limit our consideration to constraints 
on single parameters, fully marginalized over all the other parameters. First we calculate the 95% error bars with only the 
statistical errors. Then, for each systematic effect, we change our handling of the effect as described more completely below for 
each case. When we do this, the 95% confidence intervals move around somewhat. We define the upper systematic error to be 
the maximum upper limit of the confidence interval allowed by the various changes minus the nominal upper limit with only the 
statistical errors. Likewise the lower systematic error is the lower limit with only the statistical errors minus the minimum allowed 
lower limit. Finally, we (conservatively) estimate the total error as the sum of the statistical errors and each of the systematic 
errors added linearly, not in quadrature. 

For the residual PSF (B) systematic, we implement the same technique we used in iJarvis et al.\ J2003I) . namely running the 
analysis with the S-mode contamination added to and subtracted from the i?-mode signal. Most types of contamination either 
add power (roughly) equally to the E and B modes, in which case the subtraction is appropriate, or mix power between the two 
modes while conserving total signal, in which case the addition is appropriate. We allow for both possibilities to estimate how 



4 Wayne Hu, private communication 

5 It is for some of these high w a models that the ISW effect may be important, requiring a more careful analysis to determine what portion of the nominally 
allowed region is really ruled out. 
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TABLE 1 

Fully Marginalized Parameter Constraints (95% c.l.) 



Parameter 


CMB+SN 


+Lensing 


Estimated Systematic Errors 


Final 




Constraint 


(stat. only) 


B 


CAL 


Z 


NL 


Estimate 


W = — 1, W a : 


= prior: 
















301+ 007 
0.953±°;S 


0.256+™! 
0-812±g;Sg 


+0.005 
-0.001 
+0.028 
-0.006 


+0.008 
-0.003 
+0.019 
-0.006 


+0.008 
-0.005 
+0.020 
-0.007 


+0.000 
-0.000 
+0.001 
-0.001 


0.256+° ;° 7 
0.8121°;" 


-3 < w < 0, 


w = prior: 














W 


285+ - 08 
u,z ™-0 .07 

0.846+°-?? 
-0.935±*£ 


0-2541°;" 
0.790l°;iJ 
-o.894lg;lS 


+0.005 
-0.001 
+0.022 
-0.015 
+0.006 
-0.016 


+0.006 
-0.005 
+0.018 
-0.012 
+0.005 
-0.017 


+0.006 
-0.001 
+0.019 
-0.012 
+0.005 
-0.019 


+0.000 
-0.000 
+ 0.001 
-0.001 
+0.001 
-0.000 


254+ ' 07 
UJJ *-0 .05 

0.790±g-^ 
-0.894l°; 2 f 


-8 <w < 8, 


-8 < MJ a < 8, 


Qdb(z = 1100) < 0.5 


prior: 








£2 m 

^8 
W 
W a 
w(z = 0.4) 


29+ 011 
u ' za -0.06 

u - ' 9 -0.35 

1 17 +0.50 
i - i '-0.57 

i.36i;-q 

-0 90+ ' 20 


29+ ' 11 
u - zy -0.07 

74+0-13 
u -'^-0.17 

_1 14+0- 41 

1 48+0-53 
J -*°-1.57 
n S7 +0.13 


+0.009 
-0.000 
+0.044 
-0.029 
+0.008 
-0.044 
+0.101 
-0.258 
+0.005 
-0.047 


+0.001 
-0.000 
+0.013 
-0.016 
+0.006 
-0.020 
+0.021 
-0.055 
+0.000 
-0.028 


+0.009 
-0.001 
+0.024 
-0.000 
+0.016 
-0.004 
+0.038 
-0.124 
+0.002 
-0.007 


+0.000 
-0.000 
+0.000 
-0.000 
+0.002 
-0.000 
+0.013 
-0.082 
+0.000 
-0.001 


29+ 13 
u - zy -0.07 

74+0- 21 
u -'^-0.22 
_1 14 +0-45 
x - lyi -0.60 

i-«±S;S 

-0.87lg : S 



the contamination could be affecting the cosmological fits. 

The calibration (CAL) uncertainty includes errors in the dilution calculation, the responsivity formula, and possibly biases in 
the shape measurements. This systematic was the subject of significant discussion at the recent IAU symposium 225 6 . There 
seem to be calibration differences of order 5% in shear estimators between differe nt method s. Tests with simulated images with 
known shears indicate that we have calibration errors of less than 2% dHevmans et a/120 05). but we allow for ± 5% in our shear 
values as a conservative estimate of this systematic. 

For the redshift calibration (Z) of our s urvey, there are t wo public redshift surveys with depths similar to our observations : 
the Caltech Faint Gal axy Redshift Surve y (ICohen et al.l200(jt) (CRS), and the Canada-France Redshift Survey dLillv et all 19951) 
(CFRS). We argue in lJarvis et al.\ 112003) that the CRS is a better choice, since it is more complete in the R filter band pass used 
for our observations. However, switching to the CFRS distribution allows us to estimate the uncertainties due to the redshift 
calibration. Also, since we only have one other survey to use, we cannot run symmetric plus and minus versions of this test. 
So when the 95% confidence limit moved inward for a value, we take the absolute value of the change as the measure of the 
systematic error, since a different redshift survey might hav e moved the limit a similar amount outward. 

Finally, for the non-lin ear predictions (NL), we used the Smi th et ali (|2003) model, which were an improvement over that of 
iPeacock & Doddsl (n"996). Switching back to the older model should give us a rough (over-)estimate of the remaining uncertainties 
due to the non-linear modelling. Again, we cannot run symmetric tests, so we take the absolute value of any change as a measure 
of the systematic error. 

Since this technique is non-standard and may be confusing, an example with the actual values might help explain it. For 
the CAL test with the ACDM prior, when we multiplied the shear data by 1.05 (for the +5% test), the upper limit of the 95% 
confidence interval for er 8 moved from 0.910 to 0.929, and increase of 0.019. This is our estimate of the positive systematic error. 
Similarly, in the -5% CAL test, the lower limit decreased, providing the negative systematic error for er 8 . 

In Tabled we present the statistical and systematic error estimates for il m , erg, wq and w a , for each of our dark energy 
priors. The first two columns present the estimates for each value with and without the lensing data, marginalized over all 
other parameters, and quoting only the statistical errors. The next four columns show the estimated systematic error due to each 
effect listed above. The final column includes the total uncertainty with the systematic errors added linearly with the statistical 
uncertainty. Note that we make no attempt to estimate the systematic errors present in the CMB or SN data. 

It is apparent that the systematic uncertainties in our survey are smaller than the statistical uncertainties for each of the above 
cosmological parameters. However, there is still significant room for improvement in all of these sources of systematic errors. 
Future lensing surveys which expect to reduce the statistical uncertainties by a factor of order 10 will need to address these 
systematics so that they do not dominate the final error budget. And while the systematic errors due to the non-linear modelling 
are essentially completely negligible for our survey, they will be more significant for lensing surveys with smaller fields than 
ours. 

5. DISCUSSION 

We have used our measurement of the shear two-point correlations to constrain the clustering of mass at redshifts z ~ 0.3 and 
the density and equation of state of dark energy. This has been done by combining the lensing information with the CMB and 
supernovae data. The three probes are sufficiently complementary that the joint contraints are significantly better than from any 
one or two methods. 

6 July 19-23, 2004, Lausanne, Switzerland 
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We find that the primary sytematic effects on our lensing data are the redshift distribution of the galaxies, the overall calibration 
of the shear estimates, and the systematic error due to the coherent PSF anisotropy. With our new analysis, the total effect of 
these three systematics is smaller than the statistical errors. The PCA techniq ue has substantially reduced the contribution of the 
B-mode, which used to be the dominant systematic error (tJarvis et a/.ll2003l) . We are continuing to work on methods to reduce 
this and the calibration uncertainty. 

Improving the redshift distribution would require more data: either a larger redshift survey of similar depth as our data, or 
imaging the galaxies in three or four other filters to measure photometric redshifts. As discussed below, this second option would 
also allow us to bin the galaxies and make tomographic measurements. 

We have necessarily made some choices of priors and data sets in our paramet er analysis. The datasets we have used in addition 
to the lensing data a re the WMAP first year extended data (Verde et al. 2003) (which also include CBI and ACBAR data) and 
the Riess et al. (2004) Type la Supernovae data. We have assumed that the universe is spatially flat; weakening this assumption 
significantly weakens constraints on dark energy, especially if w is allowed to vary in time. We have also assumed no tensor 
contribution to the CMB power spectrum, that the primordial power spectrum is an exact power law (no running), a nd we have 
negle cted the effects of massive neutrinos on the power spectrum. Current upper limits from cosmology (see e.g. Seli ak et all 
2005) are below 1 electron volt for the the sum of neutrino masses. Allowing for massive neutrinos could lead to (at most) a few 
percent increase in our estimated er 8 , as the presence of massive neutrinos suppresses the power spectrum on scales that affect 
our observed shear correlations, but would not lead to interesting constraints on the neutrino mass. 

Combining our data with CMB and SN data, our investigation of dark energy models show no evidence for the dark energy 
being different from a standard cosmological constant. Constant w models are consistent with w = — 1, and variable w models 
are consistent with w a = 0. The constraints on dark energy from weak lensing come from a range of redshifts centered at z ~ 0.3, 
and extending by about 0.2 in redshift on either side. Thus the measurements of its density and of w should be interpreted with 
this redshift range in mind. When we combine our data with supernova data, we are using information from different redshifts, 
and the combined data have a pivot redshift of about 0.4, where w is best measured. At this redshift, we also find that w is 
consistent with —1. 

Our analysis can be compared with other recent work that combines CMB and supern ova data with galaxy clustering, the 
abundance of galaxy clusters, the clustering of the Lyman alpha forest and other pr obes ( Soer geTef a/j l2003: Well er & Lewis! 
120031 iTegmark ef"alll2004t IWang & Te gmark 2004; Seliak et al. 2005; Rapet tTef q/.ll2b 05). It is a powerful consistency check 
that these different methods appear to agree in their conclusions. It is interesting to compare the different redshi ft ranges p robed 
by the se methods, and explore constraints on the time dependence of the equation of state (e. e. lLinderl2003l iHuterer & Cooravl 

The prospects for constraining dark energy with future lensing surveys are very interesting. With a well designed survey and 
the recent analysis techniques, it can be hoped that systematic errors would stay at levels comparable to statistical errors. The 
addition of tomographic informatio n with photometric redshifts would allow for significantly better constraints on cosmological 
parameters from lensing alone dHull!999h . The use of three point correlati ons would allow for some independent checks on 
systematics as well as improved constraints on cosmological parameters (see lPen et al.ll2003t iJarvis. Bernstein & Jai n 2004 for 
detections of the lensing skewness, and Takada & Jain 2004 for forecasts). Thus even with a survey of size similar to ours, 
significant improvements in parameter measurements are possible. With future surveys that will cover a significant fraction of 
the sky, weak lensing should allow for very precise measurements of the mass power spectrum, the dark energy density and its 
evolution. 

We thank Wayne Hu, Eric Linder, Masahiro Takada and Martin White for helpful discussions. We are grateful to Licia Verde 
and the WMAP team for making available their Markov chains and to Adam Riess, John Tonry and the High-z Supernova team 
for making available their data and likelihood code. We also thank the anonymous referee for useful comments. This work is 
supported in part by NASA grant NAG5- 10924, NSF grant AST-0236702, and a Keck foundation grant. 
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